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Abstract 

We analyze the statistical properties of three-dimensional {2>d) turbulence in a rotating fluid. 
To this end we introduce a generating functional to study the statistical properties of the velocity 
field V. We obtain the master equation from the Navier-Stokes equation in a rotating frame and 
thence a set of exact hierarchical equations for the velocity structure functions for arbitrary angular 
velocity f2. In particular we obtain the differential forms for the analogs of the well-known von 
Karman-Howarth relation for Sd fluid turbulence. We examine their behavior in the limit of large 
rotation. Our results clearly suggest dissimilar statistical behavior and scaling along directions 
parallel and perpendicular to fi. The hierarchical relations yield strong evidence that the nature 
of the flows for large rotation is not identical to pure two-dimensional flows. To complement these 
results, by using an effective model in the small-0 limit, within a one-loop approximation, we 
show that the equal-time correlation of the velocity components parallel to displays Kolmogorov 
scaling q~^/^ , where as for all other components, the equal-time correlators scale as q^^ in the 
inertial range where q is a wavevector in 3(f. Our results are generally testable in experiments 
and/or direct numerical simulations of the Navier-Stokes equation in a rotating frame. 
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I. INTRODUCTION 



Near the second order phase transition equihbrium systems exhibit scahng behavior for 
thermodynamic functions and correlations. These are characterized by certain scahng ex- 
ponents which depend on the spatial dimension d and the symmetry of the order parameter 
characterizing the phase transition. These, however, do not depend on the parameters spec- 
ifying the Hamiltonian [l| . Time dependent correlation functions, characterized by dynamic 
scaling exponents also show similar universality [;2|. These standard universal properties 
of equilibrium critical dynamics are fairly robust with respect to perturbations violating 
detailed balance j^. In contrast, truly nonequilibrium systems, like fluid and magneto- 
hydrodynamic turbulence, surface growth etc., are described by appropriate equations of 
motion and exhibit much richer universal behavior. Non-equlibrium systems tend to be 
more sensitive on the parameters that appear in the equations of motion. For example, one 
finds that for the Kardar-Parisi- Zhang equation, anisotropic perturbations are relevant in 
d > 2 spatial dimensions, leading to rich phenomena that include novel universality classes 



and the possibility of first-order phase transitions and multicritical behavior 



Turbulence in fluid, described by the Navier-Stokes equation for the evolution of 

the velocity field v, is a good candidate of systems out of equilibrium due to the external 
drive acting on the system. Statistically steady fluid turbulence in three- {3d) and two- {2d) 
dimensions show markedly different behavior: In 3d, homogeneous and isotropic turbulence 
is characterized by a set of multiscaling exponents for the structure functions (see below) 
for distance r in the inertial range between the forcing scale L and the dissipation scale rj^ 
(i.e., rjd r <^ L), and forward cascade (from small to large wavenumbers) of the energy. 
Turbulence in 2c? shows an inverse cascade of kinetic energy from the energy-injection scale to 
larger length scales and a direct cascade in which the enstrophy cascades towards small length 

l; 

scales 13]; in many physical realizations of 2d turbulence, there is an air-drag- induced friction. 
In this direct-cascade regime, velocity structure functions show simple scaling but their 
vorticity counterparts exhibit multiscaling jsl], with exponents that depend on the friction. 
In a rotating fluid isotropic symmetry is broken by the global rotation. How the breakdown 
of rotational invariance affects the universal properties remains a very important theoretical 
question. These studies are also important for geophysical flows, e.g., flows in ocean and 
atmosphere. Ref. 9|, in numerical simulations of forced rotating incompressible turbulence 
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within a periodic box of small aspect ratio, showed that above a critical rotation, 3d forcing 
leads to a 2d inverse cascade. Further Ref. [lO|, in a spectral approach to rotating turbulence 
applied to a specific eddy damped quasi-normal Markovian model, showed a trend towards 
two dimensional behavior in presence of rotation. Ref. in a helical decomposition of 
the Navier-Stokes equation demonstrated similar trends of two dimensionalization of energy 
transfer. In recent studies using a shell model for rotating turbulent fluid, the authors showed 
how quasi two-dimensional behavior emerge as rotation speed f2 increases [12I. Recent Direct 
Numerical Study (DNS) [3] for rotating fluid turbulence suggest that the energy spectra 
for the velocity components perpendicular to Q, scales as gj^ in the inertial range, with 
wavevector q_|_ being perpendicular to Q. It is not fully established whether the statistical 



m 



nature of the flow for large-fi is truly two-dimensional or not. In this context, Ref. 
a helical decomposition argued that at the lowest order rotating turbulence is not the same 
as 2d turbulence. 

Since a complete description of fully developed homogeneous and isotropic 3d turbu- 
lence requires enumeration of all the multiscaling exponents, it is important to find out the 
nature of multiscaling in the presence of rotation. Until now, there is no theoretical per- 
turbative calculational framework to obtain these multiscaling exponents within controlled 
perturbative approximations. However, the multiscaling exponent for the third-order struc- 
ture function for the longitudinal component of the velocity field is well-known: Defining 
Av(2;) = v(xi + x) — v(xi), one has in the inertial range 

([At;(|r|).ff) = -ler, (1) 

where e is the total energy dissipation (see below') and f = r/r. This is the well-known von 



Karman-Howarth 4/5-law of 3d fluid turbulence ISj. The corresponding differential form is 
given by 

Vj{Av^Avj) =4e. (2) 

Subsequently, an analog of Eq. ([2]), involving mixed correlation tensor of the velocity and 
vorticity, has been obtained for helical turbulence [16j. In particular they find correlation 

{Vi{x)vj{x.)wj{^ + r)) = --^ri, (3) 

It remains to find out analogs of these relations as above in the presence of rotation, which 
introduces helicity in the system. 



Our main results in this article are (obtained by using two different strategies): 

• By using the Navier-Stokes equation in a rotating frame, we set up the master equa- 
tion for the generating functional of the probability distribution of the velocity field 
differences. We use this to obtain the exact hierarchical relations between different 
order structure functions of velocity components. In particular we obtain a set of 
differential forms for the analog of the well-known von Karman-Howarth relation for 
non-rotating 3d fluid turbulence. These relations are, however, non-integrable and 
involve contributions from pressure, unlike non-rotating 3d turbulence. Thus we find 
that, unlike the case of fluid turbulence in an inertial frame, there are no closed 
equations for the third order structure functions and hence no simple analog of the 
well-known von Karman-Howarth relation for inertial frame fluid turbulence in rotat- 
ing turbulence. A combination of these relations yield a compact differential form of 
von Karman-Howarth like relation, which although identical to its isotropic analog 
cannot be integrated owing to the underlying anisotropy. We examine the relations 
in the limit of large-f2. In this limit, we obtain simple relations connecting velocity 
structure functions with mixed structure functions, constructed out of gradients of 
pressure differences and components of velocity differences. Finally, our results here 
are illustrative of the differences between the statistical properties of 3d turbulent flows 
in the presence of large rotation and those of pure 2d turbulence. 

• From the Navier-Stokes equation in presence of rotation we define an effective model 
which is expected to be valid in the small-fi limit. We apply mode-coupling methods on 
this effective model to calculate the equal-time two-point correlation functions of the 
velocity field components Vi, i = x,y,z. We find that for i = z, {\vz{ci,t)\'^\) ~ q~^^/^, 
where as (f j(q, t)f j(— q, t)) ~ q~^, i,j ^ z, where fl is along the 2;-direction and 
q is a 3(i wavevector. Thus, correlations involving or Vy exhibit scaling ~ q^^ 
demonstrating the anisotropic nature of rotating turbulence. 

The remaining part of the paper is organized as follows: In Sec. [Illwe write down the Navier- 
Stokes equation in the presence of a global rotation. In Sec. IIIII we set up the hierarchical 
relatione between equal-time structure functions of different order and the effects of rotation 
on them. We then define our effective model to study rotating turbulence in Sec. IIV[ We 
use it to perform a one-loop self-consistent (OLSC) calculation on the effective model to 



calculate various equal-time velocity two-point correlation functions and enumerate their 
scaling in the inertial range. We examine the limit of large-Jl. Finally, we summarize and 
discuss our results in Sec. El 



II. EQUATION OF MOTION 

The Navier-Stokes equation in presence of a global rotation O = Qz is given, in the 
rotating frame, by 

— + 2(rj X v)+v.Vv = ---^ + i^oW + f (4) 
at p 

with V.v = (conclusion of incompressibility) . In Eq.(|4]), p* is the effective pressure which 
includes the centrifugal force (through an effective potential ||f2 x r|^), z/q is the kinematic 
viscosity and f is a large-scale external forcing function, required to maintain a statistical 
steady-state. The term 2(17 x v) is the Coriolis force. Note that Eq. (jl]) is invariant under 

z -z, -v^. (5) 

However, the symmetry under — t- — r_|_, v_|_ — t- — v_i_, which is present in the Navier-Stokes 
equation without rotation, is broken in Eq. (jlj) due to rotation, where r_L = {x,y), v± = 
{vx, Vy). Thus rotation breaks the invariance under in-plane parity inversion. In a nonequlib- 
rium system there is no particular relation between the noise variance and the dissipation 
coefficient, unlike in equilibrium systems where such a relation exists due to the Fluctuation 
Dissipation Theorem [l7|. A promising starting point for theoretical/analytical studies on 
homogeneous and isotropic fluid turbulence is the randomly forced Navier-Stokes model, 
where the forcing function f is a Gaussian random force whose spatial Fourier transform 



f(q, t) has zero mean and covariance 18 1 
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{hMa-ci,-u)) = p-,M)^y (6) 

in (i-dimensions where Pij = 6ij — is the transverse projection operator in the Fourier 
space, and q and u are wave vector and frequency respectively. This was used in Ref. [1 
to calculate various universal quantities associated with 3d turbulence. The Model A (with 
y = 2 — d) and Model B (with y = 4 — d) of Ref. [l9| may be considered as special cases 
of ([6]). In absence of any global rotation [Q = 0) in 3d, for y = 4 one obtains the famous 



Kolmogorov (K41) result for the energy spectrum: E{q) = Koe^/^q~^^^ for wavenumber k in 
the inertial range. Here Ko is the Kolmogorov's constant (a dimensionless universal number) 
and e is energy dissipation rate per unit mass which is the flux of the energy in the steady 
state. In a rotating fluid system there are two important dimensionless numbers, namely, 
the Rossby number Ro = V/{2QL) and the Ekman number Mq = vq/{2VLL'^) where V and 
L are the characteristic velocity and length scales, respectively. From Eq.(j4]) it is clear that 
the possible modifications of the statistical properties due to the rotation is essentially a 
non-linear effect, because without the nonlinearity, the rotation only causes the x— and 
the y— components of v to rotate. From a statistical mechanics point of view, quantities 
of interests are the correlations Cij{q,uj) = (f j(q, j(— q, — w)). For homogeneous and 
isotropic turbulence (i.e., for f2 = 0), Cij{q,ijj) ~ Pij{q)q~^^'~^^^g{ijj/q^) for wavevector q in 
the inerfal rauge. x and . are the roughness and the dynannc exponents respctively. , is a 
scaling function. Due to the Galilean invariance of Eq.(jl]) x + = 1 exactly [18]. It remains 
to be seen how a finite rotation affects the scaling obtained in the isotropic situation. 



III. HIERARCHICAL RELATIONS BETWEEN STRUCTURE FUNCTIONS 



In this Section we set up the exact hierarchical equations between the two-point structure 



functions of various orders by generalizing the generating-functional method of Ref. [20|. 
These equations are however not closed. These exact relations, nevertheless, reveal explicitly 
the differences between the rotating and the non-rotating cases. 

In order to find the hierarchy of equations for the equal-time structure functions in the 
statistical stationary state of the stochastically forced Navier-Stokes equation in a rotating 
frame for the velocity field vf 

Z(Ai, A2,ri,r2,t) = (exp(Ai ■ v(xi) A2 • v(x2)) = (Zq) 

= j r't;ir't;2p(vi,V2,Xi,X2,)f:)exp[Ai ■ v(xi) + A2 ■ v(x2)], (7) 

where V is the joint probability distribution of velocities vi(xi, t), V2(x2, t). The Eq. of 
motion for Z may be derived in a straightforward way: 

dZ d^Z dZ 



l]p + + + D, (8) 



where D = (z/o[Ai ■ V^v(xi) + A2 ■ V^v(x2)]^o) are the anomaly terms, Ip = ([Ai ■ Vp*(xi) -|- 
A2 ■ Vp*(x2)]2^o) are the pressure contributions, and // = ([Ai ■ f(xi) + A2 • f(x2)]Zo)- Note 



that the master equation ([8]) is not closed because of the anomaly terms D. This nontrivial 
point has been treated by a variety of approaches ranging from approximate techniques 



to rigorous studies [20|, |2l|. The problem arises when we look at the master equation in 
the limit of vanishing viscosity vq — )■ 0; here the anomaly terms produce a finite effect. For 
instance, the finiteness of the dissipation in the limit of vanishing viscosity (discussed above) 
is just produced by the anomaly terms D. In what follows we shall specify the effects of the 
anomaly terms in a more detailed way. Further, we can apply the Furutsu-Novikov-Donsker 
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23| formalism to calculate the random forcing terms. For Gaussian distributed, white- 
in-time random forces f, which are additive terms in Eq. (jl]), we obtain the the final form 
of the master equation for equal-time, two-point generating functions Z: 

dZ cP'Z dZ 

where K{r) is given by (v(r, t) ■ f(0,t)) (to be calculated by using the Furutsu-Novikov- 
Donsker formalism), which is related to the variance of f. 

In order to proceed further it is useful to apply the basic symmetries of the dynamical 
equations to simplify the structure of the master equation. We assume that statistically 
stationary turbulence has been produced under the dynamics of the stochastically forced 
Navier-Stokes Eq. (jl]). Stationarity implies 

dtZ = 0. (10) 

The Navier-Stokes equation without rotation (i.e., Eq.(|4]) with f2 = 0) is invariant under 
uniform translation. For 7^ 0, the presence of the centrifugal force x (fi x r), which 
is contained inside effective pressure p*, breaks the translational invariance. Since we are 
considering incompressible flows, pressure p* may be eliminated by using the condition 
V ■ V = 0. We obtain 

- vV = V ■ (^^ X v) + Vj(v ■ w)vj. (11) 

The resulting equation after eliminating p* then becomes 

dv ■ 

+ 2P:,,{Q, X v), + P,,v ■ Vv, = uoV\ + U (12) 

Evidently, Eq. (fT2!) is invariant under r — )■ r + Fq where Tq is a constant vector. Thus for 
the hierarchical relations between structure functions for an incompressible fluid, we impose 
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homogeneity on Z: 



Equivalently, if we define 



we can write 



Z = 2:(Ai, A2,xi - X2). (13) 



x+ = xi+X2, r = xi-X2, (14) 



(9 (9 (9 

-Z = + ^Z, (15) 



^ Z^^Z-^Z. (16) 



dx2i dxi+ dx. ^ 
Homogeneity now implies {d / dxiJ^)Z = 0. 

The 3d Navier Stokes Equation (j4]) is invariant under the Gahlean transformation 

X = X + Uot', t = t', V(X, t) = V'(x', t') + Uq. (17) 

Since Z = (exp(Ai ■ v(xi) + A2 ■ v(x2)), this Gahlean invariance implies 

Z = Z(Ai-A2,xi-X2). (18) 

If we now introduce variables 

Aq = Ai + A2, 

A = Ai - A2, (19) 

and use the same considerations as in Eqs. f|T5l) and f|T6l) for xi, X2, we see that the Galilean 
invariance is equivalent to demanding -^^Z = = Therefore the master equation 

can be conveniently rewritten in terms of the variables A and r to get the generating function 

Z(r, A) = (exp[A ■ Av]) = (Zq), (with Av = v(xi) - v(x2) ), (20) 

which obeys the following equation in steady-state: 

d'^Z dZ 

— + 4V,..-n = + + A (21) 

where we have used f2 = Qz. Due to the rotation (about the z-axis), the system is 
anisotropic and 3d rotational invariance is broken. However the system admits restricted 
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two-dimensional rotational invariance about the 2;- axis. This dictates that Z depends sep- 
arately on z and rj_ (r_L hes in the xy-plane). We define A = (A^, rji, 772), i-e., A^, rji, rj2 are 
the 2;, r and 6'-components of A, respectively. In terms of these variables Z may be expressed 

as 

Z = {exp[\^Av^{z, r_L) + r]iAvriz, r^) + r]2Av0{z, r±)]), (22) 

where Avz,Avr and Avg are z,r and 6 components of the velocity difference vector Av. 
The anomaly terms D may be written in terms of the energy dissipation rate per unit mass 
e as given below (see Appendix for details) 

+ 2r]lr]2{V^Vr){V^ve) + 2r]2K{ViVe){V^v,)]Zo). (23) 

In a statistical steady state, these terms are nonzero constants and are related to the energy 
input that is determined by the variances of the external stochastic forces. We are interested 
in the scaling properties of various structure functions in the inertial range (r <^ L = 1, where 
L is the system size), where the force correlator K{r) may be expanded in powers of r. To 
the leading order we assume the following Taylor expansion: K{r) = K{0) —Ar'^, where A is 
a numerical constant. We then neglect the term relative to contributions from D. We now 
proceed to calculate hierarchical relations between structure functions of different orders. It 
is seen easily that differentiation with respect to A^, rji, ri2 leads to various structure functions 
of different order: We write Sm,n,p = {{Av,)"'{Avr)"{Avg)P) = ^55|Sj^|A.=o,,7i=o,r;2=o- 
Now Z as a function of A2, ?7i, 772 as defined in ( l22l) satisfies the steady-state equation 

d^Z I d f dZ\ If d d \dZ ^^f d d \ ^ 

+ —7^ r^— + — -V2^ + 771— — + 4fi -7/1— + r]2^ Z 



dzdXz r± dr± \ drji J r± \ drji dr]2 J dri2 \ dr]2 drji j 

= Ip + If + D. (24) 

It will now be useful to consider the pressure contribution: The Navier-Stokes Equation in 
an inertial frame {Q = 0) is invariant under the parity transformation: r — )■ — r, v — )■ — v, 
where r is the 3d radius vector. This yields that there are no pressure contributions to the 
well-known von Karman-Howarth relation for the third order longitudinal velocity structure 



function in non-rotating turbulence 



25| . When there is a finite rotation (fi 7^ 0), there 



is no invariance under the parity transformation above, and as a result, certain pressure 
contributions are finite (see below) which are zero in the non-rotating case. In addition, the 



9 



effective pressure p* now explicitly depends on Q as given by Eq. (11 II) . and hence will be large 
in the large-fi limit. With this discussion in mind, let us now write down the hierarchical 
relations between structure functions of various orders: Separately apply and 
setting \z, rji, rj2 to zero yield 

d d 1 

— {{Av,f) + j—{AvrAv,) + —{AvrAv,) = 0, 
oz drj_ r_|_ 

^{Av.Avr) + -^{r^iAvrf) + -{{Ave)') = 0, 

— {AvAve) + -g^{AvrAvg) = 0. (25) 

Note the above relations between the second order structure functions are independent of 
f2. In order to calculate similar relations between different third order structure functions, 
we separately apply and and set A^, rji, r]2 to zero. We obtain 

^{{Av^f) + -^{{Av-^fAvr) + —{{Av^fAvr) = -4(e,), 
oz or± r_|_ 

-^{Av,{Avrf) + -^{{AvrY) + ^{{Ave)'Avr) + ^{{Avrf) - An{AveAvr) 
dAv* 

= -2{[^Avr])-A{er), 
d d 

— {Av,{Avef) + ^{Au,{Avef) + An{AveAvr) 
oz or± 

= -4(e,). (26) 

Here = h'o^ViVzY, Cr = t'o(Vifr)^, eg = I'o^V iVg)"^ , pressure difference Ap* = p*(xi) — 
p*(x2). Adding the three relations (l26l) we obtain 

^{{Av,)') + -^{{AvzfAv,) + —{{Av-zfAvr) 
oz or± r± 

+ ^{Av-z{AVrf) + ^{{AVrf) + —{{AVrf) 

Oz or r_|_ 

+ ^{Avz{Avef) + ^{Au,{Avef) + -{{AvefAv,) 
oz or± r_L 

= -2{[^Avr])-A{e). (27) 

Here (e) is the total energy dissipation rate per unit mass: (e) = {ez) + (cr) + (ee)- Relation 
(!27|) may be written in a closed form, by using incompressibility of the fluid, as 

V, [(lAvpAt;,) + {Ap*Av,)] = -4(6), (28) 
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where we have used that (^^-Af^) = due to the invariance of the Navier-Stokes equation 
dl]) under z -z, -v^. Further, |Avp = {Av^y + (Avr)'^ + (Ave)'^. Equation ([28]) is 

the differential analog of the well-known von Karman-Howarth relation of fluid turbulence 
without rotation [6]. This may be simplified further by using incompressibility: One obtains 
V j{\A\\'^Avj) = 0, which is same as in Ref. j26|. Although it is identical to the von Karman- 
Howarth relation ([2]), it cannot be integrated due to the anisotropy. Thus, in that sense, a 
simple (integral) analog of the von Karman-Howarth relation for rotating fluid turbulence 
does not exist. Equation (128|) is one of the main results of this article. When = 0, the 
pressure contributions in Eq. fl26|) vanish and consequently one obtains the result for the 
nonrotating case. Three more relations, analogous to ( l26ll may be obtained by separately 
applying o^^, on Eq. (El) and setting A^, r]i, r]2 to zero. We find 

^{{Av^fAvr) + -^{Av,{Avrf) + —{Av,{Avrf) + —{Av,{Avef) - AQ{Av,Avg) 
oz or± r± r± 

d d 

— {{Av.fAvg) + ^{Av.AvrAvg) + AVL{Av,Avr) 
oz or± 

— {Av^AVrAvg) + -^{{AVrfAvg) + —{{Ave)'AVr) + 4n[-{{AvgY) + {{AVrY}] 

= - ii^^M), (29) 

where e^e = t^o(V jVr)(y jVq) is non-zero in the presence of rotation, but zero when = (due 
to the invariance under the in-plane parity inversion). Similar to Eqs. (1261) . Eqs. ( l29l) do not 
decouple, have second order structure functions and pressure contributions. In particular, 
with the exception of the two-point structure function (in our notation defined above) 6*2,0,0 = 
((Af^)^), all other two-point velocity structure functions appear in the relations involving 
third-order structure functions (I27|) and (1291) . Despite the limitations of these relations, 
useful information may be obtained in the limit Q ^ 00 from relations (12 6 p and (12^ . Note 
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that, from Eq. ( ITT]) , in the hmit p* /Vt does not vanish, and hence we obtain 

{AveAv,) = l^{{^Ap*)Avr)=0, 
1 dAv* dAv* 

1 dAv* 
(A..A.,> = --(—Av,), 

-{(Av,f) + {(A^.,)^) = --L{^^A,.„). (30) 

These are exact relations between certain two-point velocity structure functions and two- 
point mixed structure function of pressure gradients and velocity differences in the limit of 
large rotation. Thus we find that in this limit (i) (AveAvr) = 0, (ii) second order structure 
functions (AvzAve), (AvzAvr) and the combination ((A^v)^) — ((Afg)^) are simply related 
to certain mixed structure functions involving pressure difference gradient and velocity dif- 
ference, (iii) second order structure functions 5*2,0,0 = ((Afz)^) has no simple relation with 
any mixed structure function. This is a possible indication of 6*2,0.0 having different scaling 
properties in the large-f2 limit, a fact we demonstrate below explicitly in a perturbative 
calculation by using an effective model. Finally, one may eliminate pressure p* from the 
exact hierarchical relations f l28l) and ( 129|) by using Eq. ( ITTl) . although the ensuing forms of 
Eqs. ( I26ll and ( l29l) are not particularly illuminating (we do not report their explicit forms 



here). Let us now consider our results as above in the context of the helical ana' 



Karman-Howarth relation as studied in Ref. 



16| . It may be noted that Ref. 



og 



of the von 



16| considers 



helical turbulence which is still isotropic, where as rotating turbulence as considered here 
not only has finite helicity injected into the system by the rotation, the system no longer re- 
mains isotropic due to the existence of a preferred direction (the rotation axis). As a result, 
Ref. jl6j has been able to obtain a closed form explicit analog of the von Karman-Howarth 
relation ([1]) for isotropic non-helical 3d turbulence. In contrast, we are able to obtain only 
a differential form of a von Karman-Howarth like relation for rotating turbulence. Apart 
from the anisotropy effects, yet another notable difference between our results and those 



in Ref. 16|] is the nonlocal nature of the relations for rotating turbulence. Thus, despite 



the presence of helicity in both the cases, our results are not exactly equivalent of those in 



Ref. jl6[ in any special limit. 

In the same way as above, analogous exact relations between certain third order velocity 
structure functions and mixed third order structure functions involving two factors of velocity 
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differences and one factor of pressure gradient difference are obtained in the large rotation 
limit. We write some of them below: 

- 2((At;,)2At;,) + ((At;,)3) = -l-{(^Ap*)AveAvr), 



-2{{AvefAv,) + 2((At;,)2At;,) 



1 dAv* dAv* 
.-[{^Av.Ave) + {^AveAv.)], 



-2{AvrAvgAv^) = l^[{^{Avey) +2{^Av,Avr)]. (31) 



dr I 



Relations (|3T|) may further be rewritten in terms of v only by eliminating pressure using 
Eq. (fTT]) . Again, as for the second order functions, third order function {{Avz)^) does not 
appear in these relations. This again is an indirect demonstration of anisotropic effect of 
rotation. Relations ( 125|) and ( 130|) together constitute a set of exact relations between different 
second order structure functions in the limit of i7 — ?■ oo. Similarly, relations f l2B]) . fl2^ and 
f l3ip together constitute a complete set (differential) relations involving different third order 
structure functions for rotating turbulence in the large-fi limit. One may further obtain for 
any odd positive integer n 



^dAp* 



(Av^r) = 



(32) 



for large rotation. Although, these hierarchical relations do not yield the scaling exponents 
directly, they serve as benchmarks on any theoretical (analytical or numerical) and exper- 
imental results on the structure functions. In order to progress further from Eqs. f l3Up or 
( 13T|) . one must make further approximations, see, e.g., Ref. |27|. Lastly, how do the relations 
f pOj) and fl3T]) compare with those for pure 2d turbulence? As shown in Ref. 2J], for pure 
2d turbulence, pressure does not contribute to the relations between second order structure 
functions or third order structure functions. More generally, the large-f2 limit of the steady- 
state master equation is not the same as that for d = 2: The master equation in the 
large- fi limit is 

as all other terms in Eq. f l24|) vanish in the large-O limit. The corresponding master equation 
for pure 2d turbulence in the steady state is given by [24 1 



d d 2d m d d 



2 1 



dr± dri2 r_|_ dri2 r_L dri2 drj^ 



V2 d 
r± dril 



Zo 



If2 + I 



p2 



Do 



(34) 
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where Z2 is the 2d analog of Eq. (l22i) . 7/2, /p2 and are the force contributions, pressure 
contributions and anomaly terms in 2d. Clearly, Eq. f lM|) . the master eqution for 3(i rotating 
turbulence, is different from Eq. fl34l) . the master equation for 2d non-rotating turbulence. 
This led us to generally conclude that the statistical properties of a turbulent rotating fluid 
in the large rotation limit is different from pure 2d turbulence. 

Having discussed anisotropic effects of rotation on the statistical properties of Sc? ho- 
mogeneous turbulence in terms of the hierarchical relations between structure functions of 
different orders, we now set out to calculate the scaling of the equal time two-point veloc- 
ity correlation functions (equivalently second order structure functions) explicitly within a 
one-loop perturbation calculation. 



IV. EFFECTS OF ROTATION: PERTURBATIVE ANALYSIS 



In this Section we use the Navier-Stokes equation (jlj) in a rotating frame, subject to 
stochastic forcings of the type This has a long history in homogeneous and isotropic fluid 



; see also Ref. [29| for discussions 



turbulence studies which are well-documented in Refs. 
on anisotropy and helicity using randomly stirred Navier-Stokes equation. Despite some 
well-known limitations and difficulties of the method it has succeeded in predicting several 
universal numbers and scaling exponents which are close to their experimentally obtained 
values. Note in the presence of rotation the Coriolis force introduces anisotropy even at 
the linear level in the Navier-Stokes equation [Eq.(|l])]. Increasing rotation velocity O (i.e., 
decreasing Ekman number M) is expected to modify the statistical behavior. In the Fourier 
space the Navier-Stokes Eq. (jlj) takes the form 



-iixi + VQk^)vi + z^Pijp(k)Sqt;j(q, a;i)wp(k - q, a; - Wi) - ikif 



(35) 



where k is a 3(i wavevector. Here, density p has been set to unity. The inverse of the bare 
propagator matrix Gq ^ is given by 



Gq ^ 



-iuj + vok'^ + 2P^yeyzx^ 



2Pxx^xzy^ 







2Pyy€yzx^ ~^ ~t~ 2Pyx^XZy^ 

'^Pzy^yzx^ 



"^Pzx^xzy^ 



-iu + u^k"^ 



(36) 
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Here ejjfc, k = x,y or z, is the totally antisymmetric tensor in 3d: e^yz = 1 etc. Different 
elements of Go are given in Sec. IVIIII The renormalized correlation matrix Cij(k,u) is 
formally given by 

Cij(k,u) = {vi(k,u)vj{-k,-uj) = Gim(k,u;)(/m(k, w)/p(-k, -w))Gpj(-k, -w), (37) 

where Gij(k,u) is the fully renormalized version of the bare propagator Gjjo(k, given by 
the matrix fl36|) . 

To find out non-linear effects of rotation and the renormalized correlation function, we 
now need to find out the fluctuation corrections to the different elements of Go in a systematic 
perturbation theory. It is evident that any straight forward perturbation theory will be much 
more algebraically complicated due to the non-isotropic nature of Go(k, w), reflected by the 
last term in the right hand side of Eq. fl35|) . than its isotropic counter part 1^. In order 
to circumvent this algebraic difficulty we use a modified equation as an effective equation to 
calculate the renormalization of the elements of Gq. We replace, in the last term of Eq. [35l 
Vp(k, w) by GQ{k,u)fp(k,u): Thus we obtain 

G'o'^^i + «^^iip(k)Sqt;j(q, uji)vp(k - q,u - Ui) - ikip* 

= fi- 2eijpQjGo(k, u)fp(k, u). (38) 

Here, Go is the bare propagator of Eq. ( l38l) : 

go(k,a;)= . ! ,2 - (39) 

Further the effective pressure may be eliminated by using the incompressibility constraint 
V ■ V = 0. We obtain 

{-iu + Uok'^)vi + i^Pijp(k)J:qVj{q, UJi)vp(k -q,u - Ui) 

= fi- 2Pirn(k)emjp^jGo(k,Uj)fpik,Uj) 

= M^,co). (40) 
The correlation of the effective noise (pi is zero-mean, Gaussian distributed with a variance 

(0i(k,w)0j(-k, -w)) 

= 2Do\k-y\Pij - 4P,>„e^„pGo(-k, -u)P,pik)\k-^\ 

-APim {kemnp^nGo (k, u)Pjp{\i)Do | /c | 

+SDQPimPjsemnpesrq^n^r | Gq (k, Cj) | ^P^^ (k) | /c | . (41) 
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This effective model, although not exact, retains two basic effects of rotation, namely, 
anisotropicity and breakdown of parity invariance. We now use this effective noise with a 
variance ( 14 ip together with the effective model in the incompressible limit f l40p to calculate 
the inertial range of the two-point equal time correlations of different velocity components. 
The calculational advantage of the effective model Eq. (155]) stems from the fact that its bare 
propagator is identical to Go, i-e., isotropic. Below we use Eq. (HOj) as a starting point for 
our analyses below. 

Formally, fluctuation corrections to the bare propagator is given by the Dyson equation 

G;^(k,a;) = G^^5^p{}^,u) - S„/,(k,a;), (42) 

where Gapi}^,^) is the renormalized (fluctuation corrected) propagator, G'o(k, w) is the bare 
propagator and ^^/^(k, u) is self-energy which arises due to the non-linear terms and contain 
fluctuation corrections. The self-energy at the one-loop level is given by (in terms of the 
bare propagator G'o(k, w)) 

dcoi d'^q 1 1 

^is ' " 



-A PimM J 2^ ^2^^^ ^2 + ^^^4 + j,^(k _ q)2 

xPnpsik - q)(0m(q, wi)0p(-q, -wi)) 
__A^ f dui d'^q P„p,(k-q 



2 '"'"^ ' J 27r {2tiY uf + u^k^ 
2Do\q\~y 



-iui + z/o(k - q)2 



[P^p(q) + 2Pp„(q)P^^fie. 



, I^'^l ^4 + 4P™a(q)Pp/3ea..e/3.5f^'-2^^]. (43) 

+ ^^q + v^q 

Expression (143|) shows that the self-energy S^^ has an f2-independent part which is iden- 
tical to the one- loop self-energy expression for the non- rotating case, an 0(f2) part and an 
0(f2^) part: The Oi^l) contributions to different elements of are (for details see Sec. lIXp 

S*j(k,u; = 0) ~ Pue,.i^^ / dqq''~' q-y-' , (44) 

where Cartesian indices i, j refer to x, y or z, such that the first order solution to the velocity 
field is given by T,iiejzi(f)j- It is clear that there are no 0{Q) corrections to S^^. Furthermore, 
the zero- frequency one-loop contribution to Sjj is real. The 0{Q^) contribution to S^^ has 
infra-red divergence. Demanding self-consistency and noting that there are no fluctuation 
corrections to the noise strength Do [isj], we find that the scale-dependent renormalized 
viscosity ^{k) ~ uk~^^^, where we have ignored effects of anisotropy. Noting that in the 
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long wavelength limit, the effective or renormalized viscosity diverges as z/fc~^/^ [18], where 
z/ is a (dimensional) numerical coefficient, and using it in place of bare viscosity uo in (jHj) 
we obtain 

E,^{k)^Pae,,ie/'M-\ (45) 

as the 0{il) correction to the self-energy. Here, = 2il/ (uk"^^^) is the renormalized or 
scale- dependent inverse Ekman number. Thus the ^-dependent corrections will appear as a 
series in M~^{k). With this scale-dependent self-energy T,ij{k) we now calculate the different 
elements of the two-point velocity correlation function matrix. While doing this we use the 
effective equation f HOj) together with the noise variance (HTj) which now depends upon the 
scale-dependent propagator Gij. As for the effective noise, since we have 

4>i{k,uj) = fi(k,u) - 2PimemzpGpsfs (46) 

with i = x,y, z, the ^-dependent part of the effective noise (pj has no contribution from 
G22(k, w). Explicitly calculating we find 



Gil 

Gl2 
Gi3 



-ioo + uk^/^ - 2P^yn 
-iu; + i/F/3)2 + 4||(]2^ 

2 I^yry 



{-iuj + vk^l^Y + A^Q? 

[{-iu + Z/P/3)2 + 4^n^]{-iuj + Z/P/3) ' 



Q 2Pxx^ 

~ (-iu; + z/A;2/3)2 + 4f|fi2^ 



G 
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-iio + uk^/^ + 2Pxy^ 
-iu; + i/p/3)2 + 4||f]2^ 



G 
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APxxPzyQ^ - 2P,,n{-iu + z/A;2/3 _ 2Pxyn) 



[(-iu; + z/P/3)2 + 4gfi2] 



G 



33 



—iu) + z/A;2/3 ' 

Gsi = = G32. (47) 

Let us now consider various correlation functions. We find 

n \ /I n m2\ (^^(k, a;)0^(-k, -w)) 
a.(k,u;) = {\v,ik,u)\ ) = ^2 ^ ^2fc4/3 • (48) 

Thus, the equal time correlator Czzik,t = 0) ~ k'^^^'^ in the long wavelength limit. In 
contrast, all other correlators scale differently in the inertial range. Let us consider one of 
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them, namely, Cxx = (bx(k, We begin by noting that 

Vx = Gn(px + Gi2(l)y + Gis^,. (49) 

Correlation Cxx may now be calculated from ( l49l) by using the effective noise variance. A 
full calculation keeping the proper anisotropicity is algebraically difficult; hence we perform 
a scaling level calculation ignoring anisotropic amplitudes. Keeping in mind that the prop- 
agators that appear in the noise variance should be renormalized propagators, we find that, 
since (pi does not receive any contribution from Gzz, we find 

CxxiKu) = {\Gu<Px + G,2<Py + Gu<Pz\'). (50) 

Noting that the effective noise (pi involves all elements of the propagator matrix except for 
Gzz, the leading scale dependence of Cxx is given by Cxx ~ to the leading order in fl"^ 
(again we have ignored the anisotropy). In a similar way, Cyy also scales as k~^. We have 
already shown above that Czz scales as k~^^^^. Thus Cxx and Cyy scales very differently 
from Czz- The latter scales same as in 3d, where as Cxx and Cyy ~ is less steep than 
Czz- For each of these we however only extract the scale-dependence and no result on 
the dependences on kz and k separately. Let us see what our results above mean for one- 
dimensional energy spectra. Since ([^^(k, t)^) scales differently from (|t>j(k, t)p), j = x,y, it 
is natural to define two one-dimensional energy spectra: E\\{k) ~ k'^{\vz(k,ty) ~ k~^/^ and 
E± ~ k'^{\vj(k,t)\'^) ~ k~^, j = x,y, where the factor k'^ appears due to the isotropic phase 
factor in 3d (since we ignored the anisotropic coefficient we have used the 3d isotropic phase 
factor). As our scaling level results do not distinguish between different directions in the 
Fourier space, the energy spectra obtained above also do not have any anisotropy. Finally, 
in a related issue, our results here provide for a naturally occurring example of a situation 
where scaling are different in different directions, see, e.g., Ref. j3o| and references therein 
for examples of models exhibiting such 
results: Recent direct numerical studies 



behavior. Let us now compare with the existing 
13| of the forced Navier-Stokes equation @ for a 
rotating fluid for large rotation yield for one-dimensional spectrum E(k±) ~ k^, where k^ is 



perpendicular to O. Similarly numerical solutions of shell model equations 12| for rotating 
turbulence yield for energy spectrum E{k) ~ k'"^ for large rotation (shell models, being 
one- dimensional models, do not distinguish between different directions). Weak turbulence 
theory [l^ approach as well suggests a spectrum ~ kj_'^. In contrast, our one- loop mode- 
coupling calculations for small- fl yields spectra k~^, less steep than those obtained in earlier 
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numerical approaches. However, two things should be kept in mind while comparing our 
results with results obtained from other approaches: (i) our perturbative calculation and 
the effective equation are valid only for small-fi, (ii) we did not distinguish between different 
directions (due to intractable and complicated nature of the underlying one- loop integrals), 
and our consequent usage of the 3d phase factor is actually not appropriate. An anisotropic 
phase factor is likely to change the scaling of energy spectra here. These make direct 
comparisons with other existing results difficult. 



V. CONCLUSION 



In this article we have analyzed the effects of rotation on the scaling properties of 3d 
homogeneous incompressible turbulence. We have used a two-pronged strategy: First of 
all, we set up hierarchical relations between structure functions of different orders by using 



an approach used in Refs. 2^, |25||. Unlike the isotropic case, there are no closed relations 



involving the third-order structure functions. Moreover, second-order structure functions 
appear in those relations. Thus simple analogs of the well-known von Karman-Howarth 
relation for isotropic fluid turbulence do not exist here. Furthermore, mixed second-order 
structure functions made of pressure gradients and differences of velocity components appear 
in these relations. We are able obtain a differential form analog of the von Karman-Howarth 
relation of non-rotating fluid turbulence. However, this is non-integrable. All these features 
are in contrast to the results for isotropic turbulence {fl = 0). In the limit of large-Jl 
these relations yield exact relations between certain second-order velocity structure functions 
and mixed second order structure functions of pressure gradient and differences of velocity 
components. The overall structures of the hierarchical relations suggest that the scaling 
properties of velocity structure functions involving one or more in-plane (plane perpendicular 
to the axis of rotation) is expected to be very different from 5*2,0,0 = {{^Vz)"^). In a similar 
way, we are able to derive exact relations between third-order velocity structure functions 
and mixed third-order structure functions involving two factors of velocity differences and 
one factor of pressure gradient difference in the limit of large rotation. Again, similar to the 
case of the second-order structure functions third-order structure function S'3,0,0 = ((^"^2)^) 
does not appear in any of these relations, reinforcing further the effects of anisotropy due 
to rotation. We then ask: Are the turbulent flows at Q ^ 00 statistically same as pure 2d 
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turbulence? As we have discussed above, the two are not identical, a fact which is clearly 
brought out by the respective structure function hierarchies in the two cases. Thus our 
conclusions are in agreement with that of Ref. [l^. With these exact relations at hand, 
we then embark upon an explicit calculation of the scaling of different two-point equal 
time velocity correlation functions. In order to simplify the ensuing calculations we use an 
effective model and use a one-loop approximation for our purposes. The resulting one-loop 
integrals are complicated due to the anisotropic nature of the system. Treating them at 
the scaling level, i.e., ignoring anisotropy, we obtain (|t'2(k, ~ k^^^^^, identical to the 
result in three dimensions. In contrast all other two-point equal time velocity correlators 
display a scaling ~ k~^. The latter result imply a three-dimensional spectra (again at the 
scaling level) ~ k~^, different from the conclusions arrived at by using other methods for 
large rotation. We would like to emphasize that our mode-coupling results are only at the 
scaling level and hold for small rotation; hence these are only illustrative of the anisotropic 
scaling due to rotation. In order to compare with numerical results or results obtained by 
other analytical means, more elaborate calculations, keeping the anisotropic coefficients of 
the one-loop integrals, should be performed. Our results, especially on the relations between 
different structure functions, may be tested by direct numerical simulations or experiments: 
One needs to calculate/measure, e.g., the two sides of the relations (126!) or (l29l and determine 
their validity. 

Our hierarchical relations are exact; however they are not closed and cannot be solved. 
Despite that they bring out two important results (albeit indirectly): (i) scaling of the two- 
point structure function 5*2,0,0 is likely to be different from the scaling of all other two-point 
structure function, (ii) statistical properties of rotating turbulent flows are not the same 
as those of pure 2d flows. In contrast, our perturbative calculations are approximate, but 
they still provide explicit results on the scaling properties of the various two-point velocity 
correlation functions which are not in contradiction with the exact hierarchical relations. 
One-loop diagrammatic calculations presented here suffer from several limitations which are 
well-documented in standard literature. Despite these difficulties we are able to obtain use- 
ful results from them and open up many new relevant and interesting questions for future 
studies. As a next step, it would be useful to account for the anisotropy in the one-loop 
integrals and find out the anisotropic scaling functions in the expressions of various corre- 
lators. Further, it would be interesting to find out perturbatively or numerically whether 
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large rotation may lead to unequal dynamic exponents for different correlation functions. 
If that happens then it would be a natural example of what is known as weak dynamical 
scaling in the literature. Until now the latter has been observed only in simple model stud- 



ies 



3l[ |. Finally some shell model studies [12] indicated that as Q increases intermittency 



corrections to the Kolmogorov's simple scaling exponents decrease and finally disappear for 
high Q. Given the fact that shell models do not distinguish between various directions, it 
would be useful to address this question by using Direct Numerical Simulations (DNS) of 
the 3d Navier Stokes equation in a rotating frame and see whether multiscaling disappears 
only for the structure functions made of in-plane velocity components and survives for the 
direction parallel to il. In general our results suggest that experiments or DNS studies 
should measure the scaling of structure functions {{Avj)"^) for j = z,r,6 with positive n. It 
should be tested whether ((Af^)") shows simple scaling for large-f2 for all n > 0. We hope 
that our results will motivate more detailed experimental work as well in the directions as 
discussed above. 
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VII. APPENDIX I: CALCULATION OF THE DISSIPATIVE ANOMALY 

We present explicit evaluation of dissipative anomaly D: We begin with 

{uoEidi exp[X^Avz + i]iAvr + r]2Avg]) 
= ^o{K['^lvz{x2) - Vlvz{xi)]to) + t'o('7i[V2^^r(x2) - Vjvr]Zo) + z/q (r/a [Vsf e (X2) - Vlveixi 

The left hand side of Eq. (1^ , upon expansion yields various structure functions of different 
order multiplied by a factor of z/q. Since structure functions themselves are finite in the limit 
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z/q — )■ 0, the left hand side of (15T!) vanishes in the hmit z/q — 0. Thus we obtain 

D = z/o(A.[V^t;,(x2) - VK(xi)]to) + z/(r?i[V^i;,(x2) - V^t'.jZo) + z/(772[V^i;e(x2) - V\ve{^i)]Z^) 

+ 2772A.z^(ViWe)(Vii;.)]Zo), (52) 
with = (ViV^Y, er = (V,^'^)^ eg = (ViVg)^. 



VIII. APPENDIX II: EFFECTIVE NOISE VARIANCE 

By definition 

{(j)i(k,uj)(j)j{-'k,-uj)) = {fi{'k,uj)fj{-'k,-uj)) 
-{fi(k, uj)2Pim(k)emnp^nGo{-k, -a;)/p(-k, -u)) 
-(2Pi„(k)e„„pil„Go(k, w)/p(k, uj)fj{-k, -w)) 

+ (2Pirrtemnpf^nGo(k, Uj) fp(k, Uj) 

x2Pjs(k)esrgnrGoi-k,u)fgi-k, -u)). (53) 
Substituting from Eq. we obtain Eq. fllTi) . 



IX. APPENDIX III 



Here we write the full expressions of the f2-independent Jq and 0{Q) part J^: 

d^q 2DoPmp{q)Onps{c[)\q\~y 



In 



-Pzmn(k) 
\2 



(27r)3 2z/oV[g' + (k-q)2] ' 
yP,^„(k) y ^4Do^^|g|-^P„p,(k-q)Pp,(q)P^^(q: 



(54) 



X 



2z/oV[g^ + (k-q)2]J ■ ^^^^ 
The above integrals are formally divergent as the external wavevector k — )■ 0. Self- 
consistency is achieved for a scale-dependent viscosity u(k)uk~^^^ and scale-dependent in- 
verse Ekman number M{k) ~ z//c^/^/2f2. 
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